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Abstract 

We present the first experimental study on the simultaneous capillary instability amongst viscous 
concentric rings suspended atop an immiscible medium. The rings ruptured upon annealing, with 
three types of phase correlation between neighboring rings. In the case of weak substrate confine¬ 
ment, the rings ruptured independently when they were sparsely distanced, but via an out-of-phase 
mode when packed closer. If the substrate confinement was strong, the rings would rupture via 
an in-phase mode, resulting in radially aligned droplets. The concentric ring geometry caused a 
competition between the phase correlation of neighboring rings and the kinetically favorable wave¬ 
length, yielding an intriguing, recursive surface pattern. This frustrated pattern formation behavior 
was accounted for by a scaling analysis. 
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Capillary instability is a commonly observed phenomenon: a slender liquid object rup¬ 
tures into a series of droplets, driven by surface/interfacial tension ( 7 ) [T|. The droplets are 
spaced at a characteristic distance, corresponding to the fastest growing wavelength (mode). 
From Tomotika’s linear stability analysis, this mode is a function of the interfacial tension 
and the cylinder-to-medium viscosity ratio [2|. When multiple cylinders are embedded in 
parallel within the same medium, the dominant mode for neighboring cylinders can become 
correlated [ 3 ]: the droplets positioned either in-phase or out-of-phase [I|. 

Instability of non-minimal shapes is fundamentally interesting. However, despite the 
rich literature on capillary instability of straight cylinders, studies on curved objects have 
been rather lacking until recent years. Pairam et al. [ 3 ] successfully created an unstable 
ring (toroid) by injecting liquid into a rotating bath of an immiscible liquid. They showed 
that the evolution of the as-formed ring was dictated by the competition between radial 
contraction and circumferential rupture [3]. Yao et al. analyzed the Stokes flow during 
the contraction [3]. Mehrabian et al. simulated both the contraction and the non-linear 
rupture of an embedded Newtonian ring pj|. The aforementioned literature suggests that 
the characteristic contraction time and rupture time predominantly scale with the medium 
viscosity and ring viscosity, respectively. Indeed, by replacing the medium with a highly 
viscoelastic material, the two timescales can be decoupled [ 8 ]. Furthermore, the stability of 
a substrate-supported liquid ring was studied both theoretically [9] and experimentally, via 
spin-coating im solvent evaporation m and pulse-laser mm, as well as ion-beam [Hi- 

Previous research has focused on a single ring; whether and how multiple closely arranged 
rings would rupture remained unclear. This could be because none of the literature methods 
were capable of creating multiple embedded rings, with well-defined dimensions and physical 
properties. In this Letter, we report the first experimental study on capillary instability 
amongst suspended concentric rings. 

The concentric rings were created by a three-step fabrication process, which we devel¬ 
oped previously [15]. Briefly, we first imprint a concentric ring pattern on a spin-coated 
poly(methyl methacrylate) (PMMA) film, via nanoimprint lithography. Next, a layer of 
polystyrene (PS) was spin-coated onto the PMMA pattern, using a selective solvent (1- 
chloropentane). The bilayer would form individual rings upon annealing. In this Letter, we 
show three representative samples (A, B and C), imprinted with different patterns. 

The cross-sectional geometry of the as-cast patterns is illustrated in Fig. [l](al), where PS 
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Figure 1. (al-a2) Schematic of the formation of suspended PS rings upon annealing. Optical image 
of concentric PS rings (w = 3.1 five i, l = 20.2 /im, Sample A) (b) as formed and (c) after rupture, 
annealed at 160 °C for the labeled durations, (d) The number of capillary waves plotted as a 
function of aspect ratio. The solid line is based on Tomotika’s theory. The dotted line is a linear 
fit to the “In-phase” data. The empty circles are adapted from reference |5j. 

mostly segregated in the PMMA trenches. Being a non-minimum shape, the pattern would 
spontaneously evolve at a temperature above the T g s of PS and PMMA. 

At first, the corrugation was leveled by the vertical Laplace pressure (P « 27r 2 yd// 2 , where 
d and l are the height and periodicity of the corrugation, respectively [Ml 113), forming PS 
rings atop PMMA. Hereafter we refer to the PS rings as the “rings”. Fig. [Tp'2) schematically 
shows the cross-section of the rings, with w, h, h r denoting the width and thickness of the 
ring, and the residual layer thickness, respectively. 

At 160 °C, the entire leveling process completed within the first minute of annealing, as 
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the flow times of both polymers under P were very short: r/pg = 802 Pa • s, and //pmma = 
1450 Pa • s, from our rheological measurements. During this process, w of the rings decreased, 
in order to balance the surface tension of PS and the interfacial tension of PS/PMMA 
(7 ps/7ps/pmma ~ 24) |l5, il8j. The R of each ring and periodicity l remained nearly constant, 
indicating negligible contraction of the rings. The only exception was the innermost ring 
in Fig. 0) evolving into a single droplet in Fig. [I]:. However, this was not due to radial 
contraction: the ring ruptured first and then rounded up (Fig. S2), which was also observed 
on fat rings by Pairam et al. [5] 

By adopting Yao & Bowick’s solution [Sj, the contraction rate for the 2 nd innermost 
ring was estimated to be on the order of 10 _3 /im/s. Before the capillary rupture time of 
~ 20 min (Fig. S2c), R would only decrease ~ 1 //in, which was negligible in comparison 
with l. The rest of the rings had even larger aspect-ratio (p = 2 R/w) and, therefore, even 
smaller contraction rate. Here wj 2 is considered equivalent to the tube radius (a) of a toroid. 
In following discussions, we will only focus on the rupture behaviors of the rings, after the 
initial leveling process. 

Fig. [lj) shows Sample A at 5 min of annealing, forming concentric rings with w ~ 3.1 /./in, 
h ~ 1.0 /rm, h r ~ 1.1 pm and l ~ 20.2 pm. All the rings had identical w , which guaranteed 
p oc R. The cross-section of each ring, represented by a ratio of w/h ~ 3, is consistent with 
that of straight filaments after the fast leveling process [El uni ED]. At this time, the rings in 
the PMMA trenches remained continuous, albeit periodic capillary fluctuations were already 
discernible (Fig. 0 3)- After 40 min, all rings had ruptured into discrete droplets (Fig. [l]:). 
Note that the ruptured segments quickly equilibrated into droplets and were kinetically 
immobilized, because collision-based coalescence rate was extremely slow. 

We plot the number of waves ( N ) from each ring as a function of p (solid squares in 
Fig. p). The relationship agrees well with Tomotika’s theory for a cylinder (with a radius 
a = w/2) embedded in an infinite medium (solid line in Fig. W, with a slope of 0.582) [21 IS] - 
Since p oc R, the linearity N oc p implies that the breakup wavelength A = 2n R/N oc R/p 
was a constant for all the rings with sufficiently large p. 

Besides A, we were also interested in the phase correlation between the rupture of neigh¬ 
boring rings. To unambiguously identify the phase correlation, we statistical analyzed the 
coordinates of all the droplets formed.We define the phase shift (</>) locally for every droplet, 
as shown in Fig. [2}t. For an arbitrary droplet X, we find its closest pair of droplets on the 
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Figure 2. (a) Schematic for defining phase shift (0); (b) Distribution of 0. The dashed line is the 
mean average of the bar heights. 

outer ring, Y and Z, and calculate central angles a and (3. We define 0 with an angular 
relationship: 0 = 2na/(3. 

For in-phase correlation, X aligns with either Y or Z along the radus (in-phase), so 0 
becomes 0 or 2i r, respectively. For out-of-phase, a = [3/2 and, therefore, 0 = ir. Fig. [2]f 
is the distribution across the entire sample, suggesting a uniform distribution: the rings 
ruptured independently. This is not surprising, given that neighboring rings are sparsely 
distanced (21/w ~ 13). Knops et al. showed that for a viscosity ratio of 0.04, the flow 
induced by capillary rupture of a cylinder extended up to ~ 10 times its radius [1]. 

In order to enhance the hydrodynamic interference between neighboring rings, we fabri¬ 
cated Sample B with reduced l, via the same procedure but with a different mold (12 /j,m 
periodicity, and a line-and-space ratio of 1). Upon annealing at 170 °C for 5 min, the surface 
leveling process was completed, resulting in a set of denser packed rings (Fig. [fijr) than Sam¬ 
ple A (Fig. |TJd) . For Sample B: w ~ 4.6 /mi, h ~ 1.4 gm, h r ~ 1.7/im and l ~ 12.1 /im. The 
cross-section w/h ~ 2.7, again, reflected the balance between the 7ps and 7ps/pmma- Most 
critically, the 2 l/w ratio for Sample B was ~ 5.3. Upon further annealing, the rings started 
to undulate and rupture (Fig. [3]b) . After 60 min, all rings had ruptured into discrete droplets 
(Fig#), whose sizes and positions remained unchanged even after 540 min. (Fig. S3) 

Similar as Sample A, the N ~ p relationship (empty triangles in Fig. 01 ) matched 
Tomotika’s theory Ei. The only exceptions are the two innermost rings, where in-plane 
relaxation dominated and reduced the number of droplets (arrows in Fig. [3p and Fig. [3};.) 
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Figure 3. Optical images of Sample B annealed at 170 °C: (a) after leveling (w = 4.6 / /m. I = 
12.1 pm) and (b, c) rupture, (d) Distribution of (j>. The dotted line is a fit to a truncated normal 
distribution within domain [0, 27 t). 


In stark contrast, Fig. [3]d suggested a unimodal distribution, peaked near tt. This unam¬ 
biguously shows that the most probable phase correlation is out-of-phase. For this sample, 
the neighboring rings were sufficiently near (21/w ~ 5.3) to interfere with each other. Based 
on recent numerical work, for a cylinder/medium viscosity ratio of ~ 1, out-of-phase corre¬ 
lation is expected for a 21/w ratio of 3 - 10) [21]. Sample A and B had a 21/w ratio of 12 
and 5.3, consistent with their non-correlated or out-of-phase mode, respectively. 

From the previous studies on straight fitments [hi si m m , the out-of-phase mode is 
the result of synchronized flow amongst neighbors: An alternation of necking and expanding 
occurred along the orthogonal direction. Therefore, N is constrained to be identical between 
neighbors. If this is also true for concentric rings , it would contradict the observed N oc p 
(Fig. Eh). We owe the observed linearity to the locality of the out-of-phase breakup, since 
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there was no indication of long-range correlation/influence across Sample B surface (Fig. [3 Jd 
and c). The correlation became more evident starting from the 6 th ring (Fig. S4). This also 
resulted in a broadened distribution of 0. 

Further decreasing the spacing between neighbors could transition the correlation into 
“in-phase”, when the axial flow started to couple amongst neighbors ©EH- This was also 
observed in sheared polymer blends |22| . From the recent numerical work EH, we expect the 
threshold of 2 l/w for out-of-phase to in-phase transition to be ~ 3, for our system (viscosity 
ratio of ~ 0.55). However, fabricating so densely-packed rings turned out rather challenging: 
Simply increasing the cast volume of PS (higher concentration or low spin speed) would only 
result in a thick top layer, which levels into a planar bilayer during annealing, as opposed 
to forming concentric rings [15]. 

We recently discovered that strongly confined straight filaments (e.g. small h r ) always 
break up in-phase, regardless of the viscosity ratio or the substrate wettability [23] . Herein, 
we fabricated substrate confined rings (Sample C). The degree of confinement can be defined 
as H/h , where H is the overall film thickness (Fig. [TJi 2 ). The smaller H/h is, the stronger 
substrate confinement is. The H/h for Sample A (Fig. 0 and B (Fig. [3]) were 2.1 and 2.2, 
respectively; both were larger than the bulk-to-confinement threshold of 2.0 [21]. Therefore, 
both cases can be considered as weakly confined. 

Fig. [4] shows Sample C (see Fig. S5 for more snapshots). After the initial surface leveling 
within 5 min (and stable up to 180 min), w ~ 8.5 /um, h ~ 1.6/xm, h r ~ 0.5 /um, l ~ 25.0/xm. 
For Sample C, H/h = 1.2, meaning substrate exerted strong confinement on the rings. 
Its w/h ~ 5.3, indicating a flattened ribbon shape, that deviated significantly from the 
equilibrium shape {w/h ~ 3) of a weakly confined thread. 

These confined rings were much more kinetically stable: they started to rupture between 
400-600 min, which was more than one order of magnitude slower than Sample A and B 
(Fig-i- The difference cannot be adequately explained by their difference in w. This is 
consistent with literature showingsuppressed capillary instability under confinement [221 ESI 
[26] . Recent work by Alvine et al. also showed that the capillary fluctuations of polymer 
melt were dramatically hindered atop a topographic Si grating EH 

Despite the slow kinetics, the rings eventually ruptured (Fig. [4]b) . However, different 
from straight filament arrays [15, [23], these droplets radially lined up. We plot the N ~ p 
scaling (diamond symbols) in Fig. The linearity again indicates constant A for all rings. 
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Figure 4. (a, b) Optical images of Sample C annealed at 170 °C for the labeled durations. Before 
undulation, w = 8.5 //,rn, l = 25.0 //in. (c) Close-up view of the undulation, (d) Schematic for the 
recursively “inserted” waves, (e) Distribution of (f>. (f) Scaling between the order of generation n 
and R. Empty diamonds represent each new wave. The dotted line is a linear fit. 

A linear fit (dotted line) shows that the slope (0.25) is less than half of Tomotika’s theory. 
This directly translates to larger wavelength and droplet size by volume conservation. 

As previously discussed, for concentric rings, N oc p (also observed in two other strongly 
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confined patterns, Fig. S6). Apparently, more waves had been “inserted” into the outer 
rings. Here we attempt to shed fight upon this process. Fig. [4|c) shows the undulation. 
The primary correlated directions are marked “1 st (generation)”, extending radially from the 
center and perpendicular to the tangential of the rings. Moving away from the center, 
more waves were “inserted” in between the primary directions. Although the PS segments 
enveloped between “1 st ” directions all had the same central angle, their arc length (also A) 
increased with R. The increased A required a gradually less favorable undulation mode and 
built up the level of frustration. When this frustration grew sufficiently large, it could be 
released by inserting an additional wave in between (marked “2 nd ” in Fig. [i]:). Similarly, the 
3 rd generation can be found at an even larger R. 

Therefore, the most energetically favorable (least amount of frustration) mode should 
correspond to the smallest A. We denote this characteristic wavelength with A c . The upper 
bound of A should be the most frustrated wavelength 2A C (on the verge of splitting up into 
two waves). Thus we have A c < A < 2A C . A c can be directly measured by identifying the 
smallest wavelength, as labeled in Fig. [4](c) . We obtained that A c = 77.8 ±11.0 pm and the 
average A = 108.9 ± 7.0 pm, which was consistent with the lower and upper bound limit. 
Tomotika’s theory provides a wavelength estimate of 45.9 pm. However, the lower bound A c 
was larger than the prediction, due to confinement-induced wavelength increase PI (26]. 

We develop a scaling relationship to capture the recursive nature of the “insertion” be¬ 
havior. As shown in Fig. [4jd, in between the n th generation envelope (radius R ni central 
angle a n , wavelength A c ), a new wave is inserted in the middle but at a larger radius R n +i- 
Therefore, VnGN: 

o n = R n a n = X c => R n = • 2 n oc 2 n , 

2 n 1 2«i 

or equivalently In R n oc n. We statistically verify the scaling against the experimentally 
observed order of generation n and In R n (Fig. [4J and Supplementary Fig. S7). The dotted 
line is a fit to the mean value of In R n for each generation (n). The excellent linearity proves 
that, despite the randomness of surface capillary instability, our simple scaling analysis was 
capable of capturing the essential “recursive” behavior of the concentric rings. 

In summary, we developed a novel procedure that allowed us to examine capillary breakup 
of concentrically arranged PS rings, suspended atop a layer of PMMA. When the substrate 
confinement was weak, the rings broke up independently if they were far apart, but via an 
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out-of-phase mode if they were sufficiently close. For both cases, the breakup wavelength 
agreed well with the prediction by Tomotika’s linear stability theory for a fully embedded 
cylinder (approximating the ring half-width as the cylinder radius). Under significant con¬ 
finement of the substrate, the rings tended to breakup via an “in-phase” mode along the 
radial direction. The unique concentric ring geometry induced strong geometric frustration, 
which yielded a self-similar morphology that could be accounted for by our scaling analysis. 
Geometric frustration associated with curvature is a fundamentally important topic. Our 
experiments can serve as a basis for correlated capillary instability among curved objects, 
which can be a powerful tool for creating unique surface patterns. 
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